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Abstract Fifty years after Ed Salpeter's seminal paper, tremendous progress both 
on the observational and theoretical sides allow a fairly accurate deter- 
mination of the Galactic IMF not only down to the hydrogen-burning 
limit but into the brown dwarf domain. The present review includes 
the most recent observations of low-mass stars and brown dwarfs to 
determine this IMF and the related Galactic mass budget. The IMF 
definitely exhibits a similar behaviour in various environments, disk, 
young and globular clusters, spheroid. Small scale dissipation of large 
scale compressible MHD turbulence seems to be the underlying trigger- 
ing mechanism for star formation. Modern simulations of compressible 
MHD turbulence yield an IMF consistent with the one derived from 
observations. 



1. Introduction 

The determination of the stellar initial mass function (IMF) is one 
of the holly grails of astrophysics. The IMF determines the baryonic 
content, the chemical enrichment and the evolution of galaxies, and thus 
the universe's light and baryonic matter evolution. The IMF provides 
also an essential diagnostic to understand the formation of stellar and 
substellar objects. In this review, we outline the current determinations 
of the IMF in different galactic environments, measuring the progress 
accomplished since Salpeter (1995) seminal paper 50 years ago. We also 
examine this IMF in the context of modern theories of star formation. A 
more complete review can be found in Chabrier (2003a) but very recent 
results are included in the present paper. 



Mass-magnitude relations 



Apart from binaries of which the mass can be determined eventually 
by use of Kepler's third law, the determination of the MF relies on the 



2 



transformation of the observed luminosity function (LF), $ = dN/dM, 
i.e. the number of stars N per absolute magnitude interval dM. This 
transformation involves the derivative of a mass-luminosity relationship, 
for a given age r, or preferentially of a mass- magnitude relationship 
(MMR), j^im\ = {^^) X which applies directly in 

the observed magnitude Mx and avoids the use of often ill-determined 
bolometric and Tcfj-color corrections. 

Figure 1 displays the comparison of the Andersen (1991) and Segransan 
et al. (2003) data in the V band with different theoretical MMRs, namely 
the paramctrizations of Kroupa, Tout & Gilmore (1993) (KTG), Reid et 
al. (2002) for My < 9 complemented by Dclfossc et al. (2000) above this 
limit and the models of Baraffe et al. (1998) (BCAH) for two isochrones. 
The KTG MMR gives an excellent parametrization of the data over the 
entire sample but fails to reproduce the flattening of the MMR near the 
low-mass end, which arises from the onset of degeneracy near the bottom 
of the main sequence (MS), yielding too steep a slope. The Delfosse et 
al. (2000) parametrization, by construction, reproduces the data in the 
My=9-17 range. For My < 9, however, the parametrization of Reid et 
al. (2002) misses a few data, but more importantly does not yield the 
correct magnitude of the Sun for its age. The BCAH models give an 
excellent representation for 0.4 Mq. Age effects due to stellar evolu- 
tion start playing a role above m ~ 0.8 M©, where the bulk of the data is 
best reproduced for an age 1 Gyr, which is consistent with a stellar pop- 
ulation belonging to the young disk {h < 100 pc). Below ra ^ 0.4 M©, 
the BCAH MMR clearly differs from the Delfosse et al. (2000) one. 
Since we know that the BCAH models overestimate the flux in the V- 
band, due to still incomplete molecular opacities, we use the Delfosse 
et al. (2000) parametrization in this domain. The difference yields a 
maximum ^ 16% discrepancy in the mass determination near My ~ 13. 
Overall, the general agreement can be considered as very good, and the 
inferred error in the derived MF is smaller than the observational error 
bars in the LF. The striking result is the amazing agreement between the 
theoretical predictions and the data in the K-band (Figure 2), a more 
appropriate band for low-mass star (LMS) detections. 

3. The disk and young cluster mass function 

3.1 The disk mass function 

A V-band nearby LF $y can be derived by combining Hipparcos 
parallax data (ESA 1997), which are essentially complete for My < 
12 at r'comp=10 pc, and the sample of nearby stars with ground-based 
parallaxes for My > 12 with rcomp=5.2 pc (Dahn et al. 1986). Such 




a sample has been reconsidered recently by Reid et al. (2002), and 
extended to r=8 pc by Reid et al. (2004). The revised sample agrees 
within ~lcj with the previous one, except for My ^ 14, where the 2 LFs 
differ at the ~2(T limit. The 5-pc LF was also obtained in the H and K 
bands by Henry & McCarthy (1990). 

The IMFs, C(log m) = ^f^^ -, derived from the <I>v and (^k LFs are 
portrayed in Figure 3 below 1 M©. Superposed to the determinations 
is the following analytical parametrization (in (log Mq)~^ pc~'^): 



{floei: m — loe 0.2)^ ~i 
~ 2 X (0 55)2 I' m<lM0 

= 0.041 m-i-35±0-3 ,m>lM0 (1) 

This IMF differs a bit from the one derived in Chabrier (2003a) since 
it is based on the revised 8-pc The difference at the low-mass end 
between the two parametrizations reflects the present uncertainty at the 
faint end of the disk LF, near the H-burning limit (spectral types ^ 
M5). Note that the field IMF is also representative of the bulge IMF 
(triangles), derived from the LF of Zoccali et al. (2000). 
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Figure 3. Disk IMF for individual Figure 4- Disk IMF for unresolved 

objects. systems. 

A fundamental advantage of the nearby LF is the identification of 
stellar companions. It is thus possible to merge the resolved objects into 
multiple systems, to calculate the magnitude of the systems and then 
derive the system IMF. The following parametrization, slightly differ- 
ent from the one derived in Chabrier (2003a) is displayed in Figure 2 
(normalyzed as eqn.(l) at 1 Mq, where all systems are resolved): 

«log m) = 0.076 X exp{- "°«";;^'g,°-^^'' }. » < 1 Me (2) 

As shown by Chabrier (2003b) and seen on the figure, this system 
IMF is in excellent agreement with the MF derived from the revised 
HST phometric LF (Zheng et al. 2001), showing that the discrepancy 
between the MF derived from the nearby LF and the one derived from 
the HST stemmed primarily from unresolved companions in the HST 
field of view. 

The brown dwarf regime. 

Many brown dwarfs (BD) have now been identified down to a few 
jupiter masses in the Galactic field with the DENIS, 2MASS and SDSS 
surveys. Since by definition BDs never reach thermal equilibrium and 
keep fading with time, comparison of the predicted BD LFs, based on 
a given IMF, with observations requires to take into account time (i.e. 
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Figure 5. Comparison of the ob- 
served and theoretical LMS, L-dwarf 
(■) and T-dwarf (a) distributions. 



Figure 6. IMF for young clusters. 
Dash- line: field system IMF (2). Ver- 
tical dotted line: H-burning limit. 



formation rate) and mass (i.e. IMF) probability distributions. In the 
present review, we proceed slightly differently from the calculations of 
Chabrier (2003a). We start from the system IMF (2) and we include a 
probability distribution for the binary frequency which decreases with 
mass. Indeed, various surveys now show that the binary fraction Xbj„ 
(and orbital separation) decreases with mass, varying from 60% 
for G and K-stars to f» 40% for early (M0-M4) M-dwarfs, to f» 20% for 
later M and L-dwarfs, correcting for undetected short period binaries, 
to ~ 10% for T-dwarfs. 

Figure 5 displays the calculated BD density distributions as a function 
of Teff, L, Mk and Mj, based on the BD cooling models developed in 
the Lyon group, and the most recent estimated LMS and BD densities 
(Gizis et al. 2000, Burgasser 2001, Cruz 2004^). The dash-hne dis- 
plays the distributions obtained with the system IMF (2) while the solid 
line corresponds to the IMF (1) for individual objects. As mentioned 
above, the distributions obtained with this latter IMF are consistent 
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with a binary frequency decreasing from ~ 50% to ~ 20%. The agree- 
ment between the theoretical calculations and the observations is very 
satisfactory, keeping in mind the remaining uncertainties in BD cooling 
theory and in accurate determinations of the observed BD T^s, Mf,oi and 
number densities. The predicted dip around Mj ~ Mk ~ 13, M^q; ~ 15 
(Chabrier (2003a)) is confirmed by the recent L-dwarf observations. 

3.2 The young cluster mass function 

Figure 6 displays the MFs derived from the observed LFs of several 
young clusters, with ages ranging from ~ 5 Myr to ~ 150 Myr, down 
to the substellar domain (sec references in Chabrier (2003a)). Note that 
some of the faintest objects in a-Or have been shown recently being field 
star contamination (McGovern et al. 2004, Burgasser et al. 2004). We 
used the MMRs from BCAH for the appropriate age in the appropriate 
observational filters. Accuracy of the BCAH models for young clusters 
has been examined carefully by Luhman et al. (2003). These obser- 
vational surveys do not resolve multiple systems, so the derived MFs 
refiect the system MFs. Superposed to these IMFs is the field system 
IMF (2). Figure 6 clearly points to a similar underlying IMF between 
young clusters and the Galactic field, except for the significantly less 
dense Taurus cluster. 

4. The globular cluster and spheroid mass 
function 

Globular clusters provide a particularly interesting test-bed to inves- 
tigate the stellar MF. They provide a homogeneous sample of MS stars 
with the same age, chemical composition and reddening, their distance 
is relatively well determined, allowing straightforward determinations 
of the stellar LF. Prom the theoretical point of view, the Baraffe et 
al. (1997) evolutionary models accurately reproduce the observed color- 
magnitude diagrams of various clusters with mctallicity [M/H] < —1.0 
both in optical and infrared colors, down to the bottom of the main 
sequence, with the limitations in the optical mentioned in §2 for more 
metal-rich clusters. As mentioned above, however, the consequences 
of this shortcoming on the determination of the MF remain modest. 
The IMFs derived by Chabrier (2003a) for several globular clusters, 
from the LFs observed with the HST, corrected for dynamical evolu- 
tion, by Paresce & DeMarchi (2000) are displayed in Fig. 7. Superposed 
to the derived IMFs is the spheroid IMF given by eqn.(20) of Chabrier 
(2003a) (short-dash line), with the characteristic mass shifted by 1 to 2 
(t's (jric = 0.22 to 0.32 M©), similar to the IMF derived by Paresce & 
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DeMarchi (2000). The disk IMF (2) is also superposed for comparison. 
We note the relative similarity of the these two IMFs, but a defficiency 
of very-low-mass objects, including BDs, in the cluster IMFs. 




Figure 7. IMF for several globular 
clusters. Dot-dash line: disk system 
IMF (2) 



Figure 8. Mass-to-light ratio for sin- 
gle stellar population in optical and 
NIR bands. 



5. Galactic implications : mass budget and mass 
to light ratio 

Integrating the IMF (1) yields the stellar and brown dwarf number- 
and mass-densities given in Table I, and the relative contributions TV = 
N(Am)/Ntot and A4 = M{Am)/Mtot, where and M denote respec- 
tively the number and mass of objects in the mass range Am. Adding up 
remnant densities (Chabrier (2003a)) yields the stellar -l-substellar con- 
tributions to the Galactic mass budget. These new determinations give 
a BD-to-star number ratio of UBD/n-^ ~ 1/3. 

Figure 8 compares the M/L ratio obtained with the present IMF (1) 
and with the Salpeter IMF (see also Bruzual & Chariot 2003). The 
present IMF yields M/L ratios in good agreement with observations 
and with expectations from CDM hierarchical simulations of spiral disk 
galaxies or from dynamical arguments (e.g. from Tully-Ficher relation), 
whereas a Salpeter IMF yields too large values (Portinari et al. 2004). 
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Table 1. Present day stellar and brown dwarf Galactic budget"'*'. 



Parameter 




Disk 


Spheroid 




Dark halo 


riBD 




2.6 X 10"^ 


3.5 X 10"^ 






Pbd 




1.0 X 10"^ 


^2.3 X 10"® 






n* 




(9.3 ± 2) X 10"^ 


< (2.4 ±0.1) X 10" 


-4 




P* 




(3.4 ±0.3) X 10"^ 


< (6.6 ±0.7) X 10" 


-5 


< 10-5 


rirem 




(0.7 ±0.1) X 10"^ 


< (2.7 ± 1.2) X 10" 


-5 




Prem 




(0.6 ±0.1) X 10"^ 


< (1.8 ±0.8) X 10" 


-5 


< lo--* 


ntot 




0.13 ±0.03 


< 3.0 X 10"^ 






Ptot 




(4.1 ±0.3) X 10"^ 


< (9.4 ± 1.0) X 10" 


-5 


< 10-* 


BD: 


M; M 


0.20; 0.02 


0.10; 0.03 






LMS(< IMq): 


N; M 


0.71; 0.68 


0.80; 0.77 






IMS(1-9M0): 


N; M 


0.03; 0.15 


0.; 0. 






WD+NS: 


JV; M 


0.06; 0.15 


0.10; 0.20 







"The number densities n are in [pc the mass densities p are in [Mq pc ^]. 



6. Star formation theory 

Although we are still far from a general paradigm for star forma- 
tion, some general properties can be considered as robust : (1) star 
formation extends well below the H-burning limit, (2) the shape of the 
IMF seems to be very similar in very diverse environments, (3) star 
formation is a rapid process, comparable to the dynamical timescale 
Tdyn = (37r/32Gp)^/^ « 1-5x10^ yr for typical star- forming molecular 
clouds, (4) the stellar IMF seems to be reminiscent of the prestellar core 
mass spectrum, suggesting that the IMF is already determined by the 
cloud clump mass distribution (sec Andre and Myers, this conference), 
(5) on large scales, the spectral line widths in molecular clouds indicate 
supersonic, super-Alfvenic conditions. All these observations point to 
a common driving mechanism for star formation, namely turbulence- 
driven fragmentation. FigTirc 10 compares the mass spectrum obtained 
from MHD simulations done recently by Padoan & Nordlund (2002) and 
Li et al. (2004) with the system IMF (2) (indeed, the prestellar cores 
correspond to multiple systems, which eventually will fragment further 
into individual objects). Although such comparisons must be considered 
with due caution before drawing any conclusion, the agreement between 
the simulations and the IMF representative of the field is amazing. Note 
that such hydrodynamical simulations form BDs in adequat numbers 
from the same fragmentation mechanism as for star formation. Vari- 
ous observations of disk accreting BDs indeed show that BDs and stars 
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form from the same underlying mechanism. Motivation for invoking 
the formation of BDs by ejection thus no longer seems to be necessary. 
Moreover, the recent observation of a wide binary BD (Luhman 2004) 
definitely contradicts the predictions of such a scenario as a dominant 
formation mechanism for BDs. This new picture thus combines turbu- 
lence, as the initial driving mechanism for fragmentation, and gravity, 
providing a natural explanation for a (scale free) power-law IMF above 
a critical mass, namely the mean thermal Jeans mass {mj), and a log- 
normal distribution below, due to the fact that only the densest cores, 
exceeding the local Jeans mass, will collapse into bound objects. Note 
that in these simulations of supersonic turbulence, only a few percents 
of the total mass end up into the collapsing cores after one dynamical 
time, solving naturally the old high efficiency problem associated with 
turbulence-driven star formation. Recent similar simulations by Li & 
Nakamura (2004), on the other hand, suggest enhanced ambipolar dif- 
fusion to occur through shock compression. 




Figure 9. Comparison of the sys- 
tem IMF (2) with the one obtained by 
Padoan & Nordlund (2002) 



Figure 10. Same with the one ob- 
tained from MHD simulations by Li et 
al. (2004). 



7. Conclusion and perspective 

In this review, we have examined the most recent determinations of 
the Galactic stellar IMF. Thanks to tremendous progress both in ob- 
servational techniques and in the theory of low-mass stars and brown 
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dwarfs, the IMF can now be determined down to a few jupiter masses, 
two orders of magnitude below the ~ 0.5 hmit of accuracy of the 
Salpctcr (1955) IMF. This IMF adequately reproduces various obser- 
vational constraints, star and BD counts, binary frequencies, galactic 
mass-to-light ratios. It is well described by a Salpeter power-law above 
~ 1 M0 rolling down into a lognormal form below this mass, with a 
characteristic mass around 0.2 Mq, although more data is needed to re- 
ally nail down this issue. Only for nearly zero-metal environments in the 
early universe do we expect a significantly larger characteristic mass, due 
to the lack of efficient cooling mechanism in the cloud (see Bromm, this 
conference). The universality of the IMF in various environments, disk, 
young and globular clusters, spheroid, points to a universal triggering 
mechanism and a dominant cooling process for star formation. Small 
scale dissipation of large scale supersonic MHD turbulence provides an 
appealing solution for this mechanism. 
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